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Abstract. This is the first of a series of papers considering symmetry 
properties of quantum systems over 2D graphs or manifolds, with continu- 
ous spins, in the spirit of the Mermin- Wagner theorem [13]. In the model 
considered here (quantum rotators) the phase space of a single spin is a 
d— dimensional torus M, and spins (or particles) are attached to sites of a 
graph (r, £) satisfying a special bi-dimensionality property. The kinetic en- 
ergy part of the Hamiltonian is minus a half of the Laplace operator —A/2 
on M. We assume that the interaction potential is C 2 -smooth and invariant 
under the action of a connected Lie group G (i.e., a Euclidean space R d or 
a torus M' of dimension d' < d) on M preserving the flat Riemannian met- 
ric. A part of our approach is to give a definition (and a construction) of a 
class of infinite- volume Gibbs states for the systems under consideration (the 
class (55). This class contains the so-called limit Gibbs states, with or without 
boundary conditions. We use ideas and techniques originated from papers 
[2], [H], [I], [22] and [7], in combination with the Feynman-Kac represen- 
tation, to prove that any state lying in the class (55 (defined in the text) is 
G-invariant. An example is given where the interaction potential is singular 
and there exists a Gibbs state which is not G-invariant. 

In the next paper under the same title we establish a similar result for 
a bosonic model where particles can jump from a vertex i 6 T to one of its 
neighbors (a generalized Hubbard model). 
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spins, symmetry group, the Feynman-Kac representation, bi-dimensional 
graphs, FK-DLR states, reduced density matrix (RDM), RDM functional, 
invariance 
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1. Introduction. Existence and invariance of a limiting 
Gibbs state 

This work had been motivated, on the one hand, by a spectacular success 
on Mermin- Wagner type theorems achieved in the past for a broad class of 
two-dimensional classical and quantum systems (see the bibliography quoted 
below) and, on the other hand, by a recognised progress in experimental 
quantum physics creating and working with thin materials like graphene. 
The main dissatisfaction with published rigorous results in this area stems 
for us from the fact that a natural class of quantum models remained un- 
covered. These are systems where the Hamiltonian contains a kinetic energy 
part given by a Laplacian. A serious problem here is that the finite-volume 
Hamiltonians are unbounded operators. As a result, the construction of the 
infinite- volume dynamical group encounters difficulties (it works fine for sim- 
plified quantum spins models (like Heisenberg's) where the phase space of a 
spin is finite-dimensional). Consequently, the KMS-definition of an infinite- 
volume Gibbs state lacks substence for this class of models, apart from the 
non-interacting case. (At least this is the situation as we know it at the 
time of writing these lines.) A consistent definition of an infinite- volume 
Gibbs state is a cornerstone for the concept of a phase transition (as a non- 
uniqueness phenomenon); it is precisely this concept that makes the Mermin- 
Wagner theorem important (and elegant). 

1.1. Bi-dimensional graphs. In the present paper we focus on Mermin 
-Wagner type result for a quantum bosonic system with continuous spins, 
over a denumerable graph (r,£) (with a vertex set V and an edge set £ C 
T x T). The graph will be assumed to satisfy a specific bi-dimensional prop- 
erty generalising properties of 'regular' lattices such as a square lattice Z 2 
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or a triangular lattice z\. Cf. Eqns (1.1.1), (1.1.2) below. (Graphene 
is clearly a regular 2D lattice; however, the whole theoretical methodology 
could be examined in the context of a more general graph with a distinct 
bi-dimensionality property.) More precisely, we assume that (T,£) has the 
property that whenever edge G £, the reversed edge G T as 

well. Furthermore, (T,£) is without multiple edges and has a bounded de- 
gree. The latter means that the number of edges (j,f) with a fixed initial or 
terminal vertex is uniformly bounded: 



sup 



max (ft O'er: (j,f) G T}, 

:{/c I': (f,j)e£}) : jer 



< oo. 



(1.1.1) 



The bi-dimensionality property is expressed in the bound 

1 



< sup 



n 



ttS(j,n): j'Gr,n = l,2,... 



< oo 



(1.1.2) 



where S(j, n) denotes the set of vertices in T at graph distance n from site 
j G T (a sphere of radius n about j): 



£(.?»= 0" G T : d(j,j') = n}. 



(1.1.3) 



(The graph distance d(j,f) = d r ,f {j, f) between sites j,f G T is defined as 
the minimal length of a path on (T,£) joining j and j'.) This implies that 
the cardinality of the ball 



A0» = O"er: d(j,j")<n}. 
grows at most quadratically in n. 



1.1.4) 



1.2. The phase space and the group action. We consider the 
following model. With each site (vertex) j G V there is associated a Hilbert 
space % realized as L/2(M,v) where M is a compact Riemannian manifold; 
v stands for the induced Riemannian volume. In this paper we assume that 
M is a d— dimensional torus R d /Z d . However, parts of the argument which 
can be easily done for a general manifold are conducted without referring to 
the specific case of the torus. (The full generalization of the main results 
for a general compact Riemannian manifold will be discussed elsewhere.) 
Physically, H is the phase space of a quantum spin 'attached' to a single site 
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of the graph and M is its classical prototype. We assume that a connected 
Lie group G is given, acting on M and preserving the flat metric on M. 
Transitivity of the action is not needed, hence G is itself a torus or a Euclidean 
space of dimension d' < d. The action is generally referred to as 

(g, x)eGxM4 gx e M. (1.2.1) 

An alternative is the additive form of writing: we represent an element g G G 
with a d-dimensional vector 

9 — 9A 

where 9 G R d ' /z d ' is a vector of dimension d! and A is a d! x d matrix with 
rational coefficient of the column rank d'. The action is then written as 

(9,x) ^x + 9mod 1. (1.2.2) 

We will use both forms: the multiplicative form (1.2.1) makes formulas 
shorter whereas the additive one is more convenient in technical calculations. 

A physical example of a system of the above type is a 'frustrated' 2D 
crystal lattice. Here some 'heavy' atoms or ions are placed at the vertices of 
a graph, and each atom possesses a light bosonic particle moving according 
to standard rules of Quantum Mechanics. A more complicated model arises 
when the number of particles is not fixed, and they can 'jump' from one 
vertex to another; see [8]. 

Another example emerges from quantum gravity: cf. [9j, [10]. Here, 
a graph is random and emerges from (random) triangulations of a 1 + 1- 
dimensional space-time complex. (The paper [§] deals with classical spins 
on random triangulations; a quantum version of the model is treated in 
[ID]-) Classical models on general graphs with a variable structure have been 
treated in a recent paper [IT]. 

If A is a finite subset in T then the phase space of the quantum system 
over A is T-L\ := H® A , the Hilbert space L 2 (M A ,v A ). Here and below the 
superscripts ®A and A mean, respectively, the tensor product of copies of 
"H = L 2 (M, f) and the Cartesian products of copies of M and v, labelled by 
sites j G A. Formally, elements of "H® A are (complex) functions 

Xa = (x(j), J G A) G M A h-> 0(x A ) G C 

considered modulo a set of t> A -measure 0, with the standard norm and the 



4 



scalar product 

i0(xA)rn^ d ^)) 1 



and 



(0i,02)= / 0i(x A )0 2 (x A ) JJv(drr(j)). 



The argument x A G M A represents a classical configuration of particles in 
A. Physically, this setting leads to a bosonic nature of the models under 
consideration. 

The action of G determines unitary operators U\(g), g G G, in "H A : 

f/ A (g)0(x A ) = 0(g^x A ) where g^x, = {g^xQ), j G A}. (1.2.3) 



1.3. The Hamiltonian of the model and assupmtions about the 
potential. A standard form of the kinetic energy operator for an individual 
spin is —A/2 where A stands for the Laplacian operator in %. We also 
assume that a two-body interaction potential is given, which is described by 
a real- valued function 

{(x',f), (*",/)) H- J(d(j',j"))V(x',x"). (1.3.1) 

In the main body of the paper we assume that the (real) function (x f , x") G 
M x M !->■ V(x', x") is of class C 2 , although in one particular result, The- 
orem 1.4, we consider an 'opposite' situation of a singular potential. (In a 
forthcoming paper, we will address in detail the case of quantum models with 
non-smooth potentials.) 

More precisely, in Theorems 1.1, 1.2, 3.1-3.2 and Corollary 3.3 below we 
assume that the function V and its first and second derivatives V X V and 
V Xl Vx 2 V satisfy the uniform bounds: x',x" G M 

\V(x',x")\, \V x V(x',x")\, |V xl V X2 lV,x")| < V. (1.3.2) 

Here x, Xi and x 2 run through the arguments x',x" G M; \ \ stands for the 
absolute value of a real scalar or the norm of a real vector, and V G (0, +oo) 
is a constant. Next, the function J : r G (0, oo) >->■ J(r) > is assumed 



5 



monotonically non-increasing with r and obeying the relation J(l) — > as 
I — > oo where 



JW,j"))i(d(f,f)>i): fer 

b"er 



J(l) = sup 

Additionally, let the interaction potential be such that 



< oo. 



(1.3.3) 



J* = sup 



U"er 



< oo. 



(1.3.4) 



Next, we assume that the function V is g-invariant: 

V(x,x') = V(gx,gx'), Vx,x' G M, g G G. 



(1.3.5) 



The Hamiltonian i^A of the system over a finite set A C T acts on func- 
tions G "H® A : given x A = (x(j), j G A) G M A , 



(// A 0)(x A ) 



5> i+ ^ wj,/))^^, 



0(x A ). (1.3.6) 



Here Aj stands for the Laplace operator in variable x(j) G M. A more 
general concept is a Hamiltonian A in the external field generated by 

a (finite or infinite) configuration xp'\A — j' G T' \ A} G M T '\ A where 
r' C T is a (finite or infinite) collection of vertices. Namely, 



(#A|x r , A 0)(x A 



jeA (j,j')eAxA 



3i X j' ' i 



X 3'> 



(j'J')eAx(T'\A) 



(1-3.7) 

Summarizing, the model considered in this paper can be called a system 
of quantum rotators on a bi-dimensional graph. 

1.4. Properties of limiting Gibbs states. Throughout the paper, 
we use a number of well-known facts (properties (i)-(iv) and (a)-(c) below) 
related to operators H\ and i?A|x r ,\ A which can be extracted, e.g., from Refs 
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: r , XA give rise 



-/3^A|x r , VA 



0, [S], [IS]) ED]- (i) Under the above assumptions, operators H A and i?A|x r , XA 
are self-adjoint (on the natural domains) in %® A , bounded from below and 
have a discrete spectrum, (ii) Moreover, V (3 > 0, i?A and H\ 

to positive-definite trace-class operators exp [— pH&\ and exp 

(iii) In turn, this gives rise to Gibbs states ip\ = tpp,A and (fApi r ,y A — f/3,A\^r'\A^ 
at temperature in volume A. These are linear positive normalized func- 
tional on the C*-algebra 23a of bounded operators in space H\. 

M^) = tr HA (R A A), Vmr ,, A (A) = tr« A (i?A|x r , VA A), A e 23a? (1.4.1) 
where 

Ra = ^ L~ — with S^,a = ti nA ( exp [ - 0H A ] ) (1.4.2) 



and 

_ exp [ - /3^ A |x r , NA ] 

^A|x r /\ A - 



-/3,A|x r , XA 

with S^|s = tr^ A (exp [ - (3H A ^ ] ) . 



;i.4.3) 



<t'\a ""AV""r L <~" il l x r'\Aj 

(iv) Let 03 stand for the C*-algebra of bounded operators in Hilbert space 
U. For A C A, the representations <B A = 23® A and 23 A o = 23® A ° identify 
03 ao with the C*-sub-algebra in 53a formed by the operators of the form A <g> 
Ia\a° where Ia\a° is the unit operator in T-La\a°- Accordingly, the restriction 
ip A of state ip a to C*-algebra 03 ao is given by 

(pf (A ) = tr WA0 (Ri°A ), A e 23 a° 5 (1-4.4) 

where 

i^° = tr„ AVAO fl A . (1.4.5) 

Clearly, operators R^ are positive-definite and have tr^ AO i?^° = 1. They 
also satisfy the compatibility property: V A C A 1 C A, 



operators R A L (> , with the same properties. 



<=tr WAlVA()J Rf. (1.4.6) 

A 

^A|x r / XA 



^A!\A° A ' 

Furthermore, in a similar fashion one can define functionals lp 1 ^"- . and 

f r'\A' 
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Below we denote by A T the net of finite subsets of T ordered by inclu- 
sion. A convenient example of an increasing sequence in this net, eventually 
covering the entire T, is formed by sets A(j,n), n = 1,2,... (balls in the 
graph distance); see (1.1.4). 

We prove in this paper the following results: 

Theorem 1.1. For all given (3 £ (0, oo) and a finite A C T, operators 
R A ° form a compact sequence in the trace-norm topology in T-L\o as A /* T. 
Furthermore, given any family of (finite or infinite) sets V = T'(A) C T and 
particle configurations xr'\A, operators Rj^j,, A a ^ so f° rm a compact sequence 
in the trace-norm topology. 

Moreover, any limiting point, R A ° , for j-R^^ , XA | IS a positive definite 
operator of trace one which possesses the following invariance property: V 

gee, 

U A0 (g)- 1 R A °U A0 (g) = R A0 . < (1.4.7) 



By invoking the diagonal process, we get a family {R A °} of positive defi- 
nite operators R A ° in 7-L A o of trace one, where A runs over finite subsets of 
T, with the following properties, (a) 3 an increasing sequence of finite sets 
A nk C T such that [JA nk = T and a sequence of sets Y' nk C Y and particle 



configurations xp^ \A nfe such that for all finite set A the convergence in the 
trace-norm holds: 

R A ° = lim R A ° |5E . (1.4.8) 

(b) V finite subsets A , A 1 of T, with A C A 1 , 

R A °=tr HAlw) R Al . (1.4.9) 

(c) Such a family defines a state (p of (that is, a linear positive normalized 
functional on) the C*-algebra quasilocal observables 2$r = (®r) • Here 
*-algebra *Bp is the inductive limit lim ind *Ba„ and superscript ~ in the 

notation (^Bp) stands for the norm completion. See [2]. 

The definition of the above state ip is that V finite A C T, 

p(A)=tm A0 R A °A. 
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Reflecting its construction, we call tp a limiting Gibbs state; Theorem 1.1 as- 
serts that the set (5° of limiting Gibbs states is non-empty. A straightforward 
corollary is 

Theorem 1.2. Any limiting Gibbs state ip G (55° has the following in- 
variance property: V finite A C T any A G 03 ao and g G G, 

<p(A) = <p(U A o(g)- 1 AU A o(g)). < (1.4.10) 



The proof of Theorem 1.1 is based on the following Lemma. 

Lemma 1.1. Let p n (x,y) be a sequence of kernels defining positive- 
definite operators R n of trace class and with trace 1 in a Hilbert space 
L2(M,i>) where v{M) < oo. Suppose there exists the following limit, uniform 
in x, y G M: 

lim p n (x, y) = p(x, y), (1.4.11) 

n— >oo 

which defines a positive-definite trace-class operator R of trace 1. Then 

lim \\R n -R\\ tT = (1.4.12) 

n— >oo 

1 /2 

where \\A\\ tT = tr (AA*) . < 

Lemma 1.1 appeared for the first time in the short note [23]. For the 
reader's convenience we give a complete proof in Section 4.3. 

Remark 1.1. As usually with Mermin- Wagner type assertions, Theorem 
1.2 does not address the issue of phase transitions, viz., uniqueness of a 
limiting Gibbs state. A matter of principle here is to determine within what 
class of states (3 D 0° the invariance property still holds true. Such a 
class is introduced in the next section; it is related to the Feynman-Kac 
representation of operator exp [ — {3H\~\ . 

Throughout the paper we adopt the following notational agreement: sym- 
bol < marks the end of a statement and symbol □ the end of a proof. 
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2. The Feynman— Kac formula and DLR equations 



2.1. The Feynman— Kac (FK) representation for the partition 
function. In this section we follow the approach developed in j6]; see also 
[I]. Our first observation is that, under the above assumptions, operator 
exp [ — /3.Ha] ac ts as an integral operator in variables xa = j G A) G 

M A and y A = {y(j), j G A) G M A : 

exp [-/3# A ]</>)(x A ) = / F[t;(di/(j))^,A(xA,yA)0(yA). (2.1.1) 

The integral kernel Xg jA (x A , y A ) admits a Feynman-Kac (FK) integral rep- 
resentation 

^, A (x A ,y A ) = / P^ AiyA (duJ A )exp [-h A (uJ A )} (2.1.2) 

x A.yA 

explained below. 

In Eqn (2.1.2), W^_ AyA stands for the Cartesian product x Wx(j),y{j)' 

Next, the Cartesian factor W x( j\ y( j\ represents the space of continuous paths 
ZJj in M, of time-length (3 and with the end-points x(j) and y(J): 

UJj : r G [0, 0\ i-> cJj(r) G M : uJj( • ) continuous, 
Wi(0) = a;(j) J Wj{P) = vij), jeA. 

Correspondingly, lJ a = (uJ,-, j G A) G M^xa.ya * s a collection of continuous 
paths aJj G ^afo-Wj); ^ e ^' ^ e sa ^ ^ A * s a P a ^ n configuration 
over A. Further, P^ AJA is the product-measure on W^ XAiYA : 

Pj AiyA (d5? A )=^pf WisW (d E ; i ) (2.1.3) 

where ^(dwj) is the (non- normalised) Wiener measure on Wa:0')>i/0') 
(the Brownian bridge in M, of time-length 0, with endpoints x(j) and y(j))- 
The measure pjv.v a y\(dw.j) is defined on a standard sigma-algebra of subsets 
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of ^ x {j),y{j) generated by cylinders, and the measure Px A , yA 011 the corre- 
sponding sigma- algebra of subsets in yA . In future we do not always 
explicitely refer to the sigma-algebras where measures under consideration 
are defined: their specification follows that of the underlying spaces. 
Finally, for a path configuration ZJ A = {ojj, j G A} over A, 

h A (uJ A ) = h"'fa,*i') (2-1-4) 

(jj')eAxA 

where h?^' (uJj,uJj') represents an integral along trajectories uJj and uJj>: 

lr hJ ' {Z: r -: r ) = J(d(j,j")) J drl/(^(r),^,(r)). (2.1.5) 



It is convenient to think that h^^'(uJj,uJji) yields the 'energy of interaction' 
between trajectories UJj and UJj/, and h A (uJ A ) equals the 'full potential energy' 
of the path configuration uJ A . 

Furthermore, the trace tr% A exp [ — f3H A ~\ (the partition function in A) 
is finite and equals S^a where 




(2.1.6) 



Consequently, operator R A from (1.4.2) (often called the density matrix (DM) 
in A) is given by its integral kernel -F/3,a(x a , y A ) (the DM kernel, DMK for 
short): 

^(xa^a) = J—Kp A {x A ,y A ). (2.1.7) 

2.2. The FK representation for the RDMK in a finite volume. 

The operator i?f from (1.4.8), (1.4.9) (referred to reduced DM, briefly, 
RDM) is determined by its own integral kernel F^ A (x A o, y A o) (the RDM 
kernel, shortly, RDMK): 

Z?A°/ \ S /3,A\Ao(x A o,y A o) _ ,,-a / o i\ 

4a( x a»Ja») = ^ , x A o,y A oGM . (2.2.1) 
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Here the quantity S /3 A \ A o(x A o, y A o) in the numerator yields a 'partial' parti- 
tion function corresponding to the partial trace tr H in (1.4.6): 

jeA\A° (2.2.2) 
xKp(x A o V Z A \ A o,y A o V z A \ao), 

where symbol V means concatenation of configurations (this notation will be 
repeatedly used below). 

It is convenient to use a brief notation dx A for the product of the Rie- 
mannian volumes x v(dx(j)). We will also omit, where possible, the argu- 

ment/index /3 from the notation (viz., by writing S A instead of Sg ;A ). The 
above representations (2.1.1)— (2.1.7) allow us to associate with Gibbs state 
y? A a probability distribution // A on the set 

W A = |J W XA , XA where W XA , XA = \JW x{j)Aj) . (2.2.3) 

x A eM A jeA 

the definition of this probability distribution is provided in forthcoming para- 
graphs. 

Pictorially, W\ is the space of collections of closed trajectories (loops) in 
M issued from and returning to (coinciding) specified endpoints; each loop 
being assigned to a site j G A. Any such loop collection can be written as a 
pair (x A ,u> A ). Here oj» a = {luj, j G A} is a collection of loops r G [0, /?] H> 
Wj(r), where ^(0) = u}j(j3) = x(j); a pair (x(j),ujj) is associated with site 
j G A. We will say that (x A ,u> A ) (and uj\ when the reference to x A is clear 
from the context) is a loop configuration over A. Note the absence of the bar 
in this notation, stressing that uj a G W\ is a loop configuration as opposite 
to a general path configuration ZU A G W^ XAiyA C Wa (again associated with 
sites j G A (see Eqn (2.2.4) below)). More precisely, when appropriate, we 
will omit the bar in the notation W X)V and W^ XA ,yA f° r x = V or x a = Ya : 

W XyX = W X:X and W XAjXA = 1^ XA , XA . 

Recall, we refer to o; A as a loop configuration and aJ A as a path configuration 
in A. Next, we set: 

W= |J W x , y andW A = |J H7 XA , yA . (2.2.4) 

x,y£M x A ,y A eM A 



12 



2.3. The FK-DLR equations in a finite volume. The aforemen- 
tioned probability distribution /i A , on space W\, is absolutely continuous 
relative to the underlying product-measure u\{= i^,a), where 

dz/ A (x A ,o; A ) = x (v((kO')) x Pxovo-)^)) 

(2.3.1) 

= dx A x P XAjXA (do; A ). 

Here the measure P x ^ tX (j)(duj) is defined as a Brownian bridge on man- 
ifold M with the starting and end point x(j). Next, the Radon-Nikodym 

derivative (probability density function) p A (x A , u? A ) := A ' — ^- is of the 

dz/ A (x A ,uJ A ) 

form 

p A (x A , « A ) = — exp [ - h A (u) A )] (2.3.2) 

"A 

where functional /i A (u; A ) has been defined in (2.1.4)-(2.1.5). It is convenient 
to treat /x A as a Gibbs probability measure for a 'classical' spin system where 
'spins' are represented by loops affiliated with sites j G A. 

To shorten the notation we will omit henceforce the argument x A and sim- 
ilar arguments from symbols like p A (x A , <*> A ), di/ A (x A ,u; A ) and d/x A (x A , u> A ), 
bearing in mind that the initial/end-point configuration x A can be recon- 
structed from the loop configuration <*> A . 

Measure ji\ defines a random field over A with realizations u; A = {ojj, j G 
A} G Wa an d has the following properties (I), (II). 

(I) satisfies the DLR equation over A; cf. Eqn (2.3.4) below. (Recall, 
A C T is a finite set.) This means the following. Given A C A, let us agree 
to write u>° for the loop configuration u> A o G W\o. Consider the partially 
integrated probability density 

Pa\u°)--= d^ A o(w A \ A o)p A (w°Vw A \ A o) (2.3.3) 

where u>° V u> A \ A o stands for the concatenation of the two loop configurations 
yielding a loop configuration over the whole of A. Cf. Eqn (2.2.2). 
Then, V set A' such that A C A' C A, the density Pa°(w°) obeys 



Pa°( w °) = / dz/ A \A'(^A\A 
Jw AXA , 

A\A' / n S A'\Ao(w ,U>a\A 

><Pa (wa\aO — ^-7- v 

-a'1^a\a'J 



A \ A ' (2.3.4) 
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(2.3.6) 



Here p^ A (u>a\aO is the partially integrated density similar to (2.3.3): 

pf- A (u>a\aO := / dz/ A /(a3 A ,)p A (a3 A / V w A \a')- (2.3.5) 
Jw A , 

Further, the quantities S A /\ A o(c<; , Wa\a') and Sa'(w A \A') are given by the 
following integrals 

S A '\Ao(w°, WA\A') — y dz/A'\A«(^A'\Ao) 

W A'\AO 

x exp [ - h A '(uj° V w A '\ao|wa\A')] 

and 

Sa'( w a\aO = y dz/ A /(a; A /)exp [ - /i A '(w A ' |<^a\A')] • (2.3.7) 

Next, given loop configurations wa' = G A'} e W A i and wa\a' = 

{wj, j G A \ A'} G M A \ A ', the functional /i a '(uja'|wa\A') in (2.3.7) is deter- 
mined by 

/i a '(wa'|wa\a') = h A '(u) A/ ) + h(u A ,\ |w A \A') (2.3.8) 

where the summand h A '(cu A i) is defined as in (2.1.4) and /i(u>a'| I w a\A') is 
given by 

h A '(u> A '||u>A\A')= ^ hM'^Uj,) (2.3.9) 
0J')eA'xA\A' 

with h?j'(ojj, Uf) as in (2.1.5). 

Finally, the functional /i A (w° Vu>a'\a°|'* , a\A') figuring in (2.3.6), for a; e 
Wap, wa' G W^a' and cl»a\a' G Wava/, is defined by similar formulas. We 
say that h A (u> A i\u} A \ A /) and h A (u>° V u>a'\a°|<*>a\a') give the values of a 'po- 
tential energy' of the loop configurations uj a > and u>° V u> A i in the external 
field generated by Wa\a'- 111 this context, S A '(w A \ A /) gives the partition func- 
tion for loop configurations over the 'intermediate volume' A' in an external 
potential field generated by the boundary condition u>a\a' G Wa\a'- Sim- 
ilarly, E A /\ A o(u)°, u>a\a') can be considered as the partition function in the 
'layer' A' \ A , with an 'external' boundary condition cl> A \a' G Wa\a' and 
an 'internal' loop configuration a; G W A o (note that u) Q enters the integral 
S A '\a°(w°, Wa\a) with its energy h A °(u) )). A straightforward fact is that 



dv A o(u°)E A /\ A o(uJ , U>a\A') = Sa'(wa\A')- 

14 



In probabilistic terms, the DLR equation (2.3.4) is equivalent to the fol- 
lowing property. Consider the conditional distribution d/i A ' A ^ A (oj°\ u?a\a') 
induced by the probability measure /ia, for a loop configuration uj° over A , 
conditioned by a loop configuration cja\a' over A\ A'. It is determined by the 

o / d/i A ° |AAA '(u>° I a; a aO 

conditional probability density p\ ' A ^ A (u>° \uja\a>) '■= — - — ; 7 — A ^ A ■ 

ais A o (ur ) 

The equivalent form of the DLR property means that this density has the 
form _ 

P f |A ^ A 'K | WAW ) = 'T 1 ; '^^ . (2.3.10) 

In fact, /ia is the only measure that satisfies the equations (2.3.5), (2.3.10). 
The name DLR (Dobrushin-Lanford-Ruelle) is widely used in the classical 
statistical mechanics; see, e.g., [5]. 

(II) The measure //a determines the RDMK F^{x°,y°). Given A°cA 
and particle configurations x°,y° G M A °, the RDMK F A °(x°,y ) is defined 
by 

F A A V,y°)= J P x o iy o(dZU°)gf(aJ°). (2.3.11) 

In turn, the functional q^°(uj°) is determined by the formula analogous to 
(2.3.4): V A' such that A'cA'c A, 



qf{u°) = / dz/ A \A'(^A\A' 



AXA ' - f -o \ (2-3.12) 

A\A' / x "A'\A°(,^ ,WA\A' 

xp A (w A \aO — ^77- x 

with quantities p^ A (a>A\A')> ^a'\a° O-^ , <*>a\A') an d Sa°(^a\a°) defined as in 
Eqns (2.3.5)-(2.3.8). (The only difference is that a loop configuration uj° 
in (2.3.7) has been replaced with a more general path configuration aJ°.) 
Kernel F i A °(x°,y ) is often referred to as the reduced DM kernel, briefly 
RDMK (more precisely, the kernel of the DM in volume A, reduced to A ). 
Accordingly, the functional q^°(uj°) can be called a reduced DM functional, 
briefly RDMF, for ZJ° e W^oy). Similarly to (2.3.4), it will be convenient to 



15 



write 



qf{uj°) = / dl/ A \A'(WA\A') 

rr a° a a' (2 ' 3,13) 

Xp^ A (wa\A')?a°' A ^ A (W°|WA\A') 

where 

A°|A\A'/_0| x S A'\Ao(w ,WA\A') ,„ 01 a 

9a (w°|wa\A') = ^-7- v • (2.3.14) 

"A'(wa\a'J 

In analogy with Eqn (2.3.10), quantity q A ' A ^ A (w° | cl»a\a' ) in (2.3.14) can be 
called a conditional RDMF for a path configuration uJ° G VF x o iy o, given a 
loop configuration cl>a\a' G W / a\a / - 

Note that the RDMF qf(oJ°) from (2.3.12), (2.3.13) satisfies the invari- 
ance relation 

(tfV) = S e G, uJ° G W x o, y o (2.3.15) 

where 

ga>° = {g^-, j G A }, with (g^)(r) = g( Wj (r)), < r < /3. (2.3.16) 

Consequently, for RDMK F^°(x°,y°) (see (2.3.11)), we have that 

Ff (x°,y°) = Ff (gx ,gy°), g G G, x°,y° G M A ° (2.3.17) 

and for the RDM Rfc° 

Ra =U A0 (g)- 1 Ri°U A0 (g), g G G. (2.3.18) 

However, we will need to consider a similar RDMF 9A°x r , XA (^°) defined 
via operator exp [ — /3i^A|x r , XA ] instead of exp [ — (3H\\ . (It can be called 
a conditional RDMF, with the boundary condition xp'\A-) The invariance 
equation 

^ A (S°) = ^ A (ff5°) (2-3.19) 

(i.e., an analog of (2.3.15)) fails; this makes the statements of Theorems 1.1 
and 1.2 non-trivial. (Of course, the covariance property 

^ r , XA (^) = < g , r , XA (g-°) 
holds true but is useless for our purpose.) 
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3. The class of Gibbs states in the infinite volume 



3.1. Definition of the class (5. The aim of this section is to define 
the invariance property (2.3.17) (and consequently, property (2.3.18)) for 
functional gf°(co> ) (and related objects F^°(x°,y ) and R A °) f° r the system 
in an 'infinite volume' (i.e., over the whole graph (r, £)). That is, we want to 
prove that functional g A °(a; ), which we call inifinite- volume RDMF, obeys 

qfiuf ) = qf(gLJ°), g G G, ZJ° G W A o. (3.1.1) 

The formal definition of infinite-volume RDMF g A (oo°) {q A (w°) for short) 
related to the system over (r, S) requires additional constructions and will 
lead us to the definition of the aforementioned class of states (5; see below. 
At this point we state that the key step is to establish an asymptotical form 
of (2.3.19) for infinite- volume conditional RDMF q A ° ' r \ A '(uJ°|a;r\A') when set 
A' is 'large enough'. In essense, we will prove that, V finite set A C T, 

_A°|r\A(n)/'_ 7 - 5 0|. , \ 

gA°|r\A(n)^o |a;nA(n)) -1. V-l-V 

uniformly in: (i) group element g G G, (ii) path configuration uJ G W\o, (iii) 
an (infinite) loop configuration u>r\A(n) G Wr\A( n ) representing an infinite- 
volume external boundary condition. Here and below A(n) = A(o,n) and 
S(n) mean the ball and the sphere of radius n (cf. (1.1.4)) around a reference 
point o G T (the choice of point o will not matter). 

In fact, functional q A ° l r \ A ( n )(uJ°|aj r \A(n)) is itself defined as the limit 

/° |rVA(n) (5? 0| ) = Um g A0| A(r)VA(«) (5? 0| ) (3 A 3) 

where, for A C A(n) C A(r), the value qf^ r) ^ n) ^\ WA(r)\A(n)) has been 
determined in Eqn (2.3.14). 

At this point it is appropriate to establish some probabilistic background. 
The Borel sigma-algebra of subsets of the loop space W is denoted by 2U. 
Given a finite subset A C T, we obtain the induced sigma-algebra of subsets 
of W\ which is denoted by QUa. Similarly, for a trajectory space W the 
sigma-algebra Q2J is defined which leads to the sigma-algebra 22Ja« of subsets 
in W\o. For A' C A, the sigma-algebra QUa' is naturally identified with a 
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sub-sigma-algebra of 2Ua which is denoted by the same symbol 2Ua')- For the 
whole graph T, we can introduce the Cartesian product Wy considered as the 
countable set of loop configurations {uj, j G T}, Uj G M; as earlier, a loop 
Uj is associated with site j G T. For a finite set Ac T, the sigma-algebra 2Ua 
can again be identified with the sigma-algebra of subsets of Wr; as before, 
it is convenient to use the same notation for both. The sigma-algebra 22?r 
is defined as the smallest sigma-algebra of subsets of Wr containing 2Ua V 
finite A C T. In a similar fashion we define the sigma-algebra 2Br\A° C 2Ur 
for a given (finite) set A C T; as before, it is naturally identified with the 
sigma-algebra of subsets in W t \ao . 

Let us now define the class (3 of states of the C*-algebra 03. As before, a 
state tp of 03 is identified with a family of RDMs R A ° = R A ° where A is an 
arbitrary finite subset of T; each R A is a positive definite operator in T-L\o 
of trace one, and the compatibility relation (1.4.9) holds true. In short, for a 
state (f G <5 the RDMs R A ° are integral operators (see (3.1.4)), with integral 
kernels F A °(x°,y°) satisfying (3.1.5)— (3.1.12), where the probability measure 
/ir obeys (3.1.13)— (3.1.17). Properties (3.1.5)— (3.1.17) are direct analogs of 
the corresponding properties of the RDMKs F J A °(x°,y°) and F^^^x , y°) 
in a finite volume A. 

Passing to the formal presentation, the RDM R A ° is determined by its 
integral kernel F A °(x°,y ): 

(i? A °0)(x°) = f d y °F A V, y °)0( y o ), x° G M A °. (3.1.4) 
Jm a " 

In turn, the RDMK F A °(x ,y ) = / A °(x°,y°) is obtained via a functional 
q A °(ZJ°) = q A ° r (cJ°) referred to as an infinite-volume RDMF: 

^ A V,y°)= J P^ iy o(d^W) 

Further, the infinite- volume RDMF for a state under consideration, should 
admit a particular representation. Namely, there exists a probability measure 
/^r = AV,r on (W r , 2U r ) such that V finite set A'cT with A C A', 

<7 A V) = / d/iP A Vr\A')/ 0|rW (^Vrw) (3.1.6) 
Jw TXA , 



(3.1.5) 
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where 

"A'( w r\A'j 

Here /xp A stands for the restriction of measure /ir to the sigma-algebra 
2Ur\A'. 

Moreover, the expressions S A /\ A o(u>°, wr\A') and S A /(u>r\A') represent, as 
before, partition functions in A'\A° and A', with the corresponding boundary 
conditions: 

E A i\ A o(p°, u? r \A') = / dz/ A /\ A o(wA'\A () ) 

w (3-1.8) 

x exp [ — h A (aJ° V cl> A /\ a o |w r \A')] 



and 



£A'(u?r\A') = y dz/ A /(w A /) exp [ - /i A '(u; A /|u> r \ A /)] . ^ j g ^ 



w 



A' 



The functional h A ' is defined by formulas similar to (2.3.8), (2.3.9): 

/i A '(u; A ,|a; r \ A ,) = h A \u A i) + h(w A /| |w r \A') (3.1.10) 

and 

/i A '(aJ° V w A /\ A o|w r \A') = h A '(uj° V w A '\a«>) /g l n \ 

+h(UJ° V w A '\ao| |w r \A') 

where 

/i(wa'| |w r \A') = 5^ hhJ ' (ujiUj') (3.1.12) 
0J')eA'x(r\A') 

and similarly for h A ' (uJ° Vu>a'\a°| |^r\A')- m turn, the terms h A ' and W^' are 
as in (2.3.8), (2.3.9). It is assumed that the series in (3.1.12) is convergent 
for /i r \ A -almost all a>r\A') £ Wr\A'- 

The functionals h A (uja>), h A (u)A'\i*>r\A'), h A (cJ° V a> A '\A°), 
h A ' (uf V u> A /\ A o|u>r\A') and h A ' (uf V u>a'\a°| l w r\A') have the same meaning 
in terms of 'energies' of loop/path configurations as before. 

The measure fj,r figuring in Eqns (3.1.6) and (3.1.8) has to satisfy the 
infinite- volume DLR equations similar to (2.3.5). Namely, consider p A °(u}°), 
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the probability density function, relative to ch/ AC V), for the measure /j, a ° , 
the restriction to the sigma-algebra 2U(A°) of measure ^r' 



rfV) = ^W "° e,1/ - (3113) 

The equations for /ip are that V finite sets A and A' where A C A' C T, 

p£V)= / d^W0l£ W Vl«^0- (3-1.14) 

Here ' r ^ A (u>°|tJr\A') is the conditional probability density for condi- 
tioned by boundary condition u>r\A' £ Wr\A' : 

A°|r\A'/ oi \ -A'\A (^°,^r\A ) /oi-ik\ 

p r (w u |wr\A') = ^— p — . (3.1.15) 

"A'l^r\A'J 

Here, as in (3.1.8) 

H A /\ A o(u; ,u> r \ A o) = / dz/ A /\ A o(w A /\ A o) 

w (3.1.16) 

W A'\A° 

x exp [ — h A (u)° V CL» A /\ A o|a; r \ A / 



and, as in (3.1.9), 

S A '(u> r \A') = J dz/ A /(w A /)exp [ - /i A '(w A /|a; r \ A /)], (3.1.17) 
w A , 

with the functional h A ' defined similarly to Eqns (3.1.10)-(3.1.12). 

Remark 3.1. We do not claim (at least in this paper and its sequel [8]) 
that the properties (3.1.4)— (3.1.17) imply that the operator R A ° is an RDM 
(positive definiteness of R A remains an open question). However, when one 
can assert (on grounds of some additional information) that a given family of 
operator {R A } obeying (3.1.4)-(3.1.17) consists of positive-definite operators 
then we can speak of a state (p G &. (Viz., this is the case of limiting Gibbs 
states discussed in Theorems 1.1 and 1.2.) As it stands, our results stated in 
Section 3.2 hold true for any family of operators R A ° for which Eqns (3.1.4)- 
(3.1.17) are fulfilled. E.g., we can claim the assertion of Theorem 1.2 for 
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any linear normalized functional on 23 defined by a family {R A } satisfying 
(3.1.4)-(3.1.17). 

Some elements of the above construction have been used in the literature; 
see, e.g., [12] and references therein. 

We will refer to Eqns (3.1.6)-(3.1.12) as an FK-DLR representation (of 
the infinite- volume RDMF q A °(u )) by a given probability measure fir, as- 
suming that /ir satisfies the infinite- volume DLR equations (3.1.14). It is 
important to stress that, unlike the case of a finite A C T, the solution to 
the infinite- volume FK-DLR equations (3.1.14) over the whole graph F may 
be, in general, non- unique. However, the family of functionals q A ° , where A 
runs through the finite subsets of F is determined uniquely provided that a 
measure fir is given, satisfies (3.1.6), (3.1.7). In accordance with the above 
scheme, this gives rise to the family of RDMKs F A and - ultimately - RDMs 
R A , for finite sets A C F, obeying the above compatibility property (1.4.9). 
The corresponding state (emerging from the probability measure fir) is de- 
noted by <^r(= ¥>(A*r)); when possible, the subscript F will be omitted. Given 
G (0, oo), the class of the measures fir = Hp,r satisfying Eqns (3.1.13) is 
denoted by C5(/3), as well as the class of related states tpr- 

In Theorem 3.1 below we establish that the class (3(/3) is non-empty V 
given (3 G (0, oo). 

As was said earlier, the infinite-volume invariance property under study 
is expressed by Eqn (3.1.1): V g G G, finite set A C F, x°,y° G M A ° and 
uJ° = {UJj, j G A } G W x o iy o, the value qf(guf) = qf(uj°). Here guJ° is as 
in (2.3.16). A similar property for the density p A0 (aJ°) has the form: V finite 
set A C T and loop configuration lJ° = {ZJj, j G A } G Wa°, 

/V)=/°(gaJ°), g G G. (3.1.18) 

The invariance properties in Eqns (3.1.1) and (3.1.18) imply that, V finite 
set A C T, the infinite- volume RDMs R A ° have the property similar to 
(1.4.9): 

R A °U AO (g) = U A o(g)R A ° (3.1.19) 
which, in terms of the corresponding state if, means (1.4.10). 

3.2. Properties of class (5. The results of this paper about state 
class (5 are summarised in Theorems 3.1-3.4 below. In Theorems 3.1-3.3 we 
assume the above conditions (1.1.1), (1.1.2), (1.3.2)-(1.3.5). 
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Theorem 3.1. For all (3 G (0,oo), the sequence of Gibbs states ipA( n ) 
contains a subsequence fA(n k ) sucn that V finite A C V and A G 25 \o, we 
have: 

lim ip k{nk) (A ) = (p(A ) 

k—too 

where state <p G Consequently, class <5(/3) is non-empty. < 

Theorem 3.2. For all (3 G (0, oo) and finite A C T, any Gibbs state 
cp G satisfies properties (3.1.1) and (3.1. 18)-(3. 1.19). < 

The invariance property can be formally extended to ground states. We 
call a state Tp (of C*-algebra 23) a ground state if there exists a sequence of 
states (p n G <&(f3 n ) with /?„ — >■ oo such that Tp = w* — lim 

Corollary 3.3. Any ground state Tp (i.e., a w*-limiting point of states 
ip n G <5(/3 n )) obeys (3.1.1) and (3.1.18)-(3.1.19). < 

Of course, Corollary 3.3 does not prove existence of ground states for the 
model under consideration. 

In a future paper we will remove the smoothness condition upon the 
potential function V (see (1.3.2)), by following the methodology from [16J. 
In this paper we note that, like the classical case (cf. [7] and the bibliography 
therein), if the condition of smoothness is violated, the symmetry property 
may be destroyed. See Theorem 3.4 below. 

Theorem 3.4. Take T = Z 2 , the regular square lattice, with distance 
d(j,f) = max - \j 2 - j' 2 \] . Take M = S 1 = G where S 1 = R/z is a 
unit circle, with a standard metric p(x, x') = min \\x — x'\, 1 — \x — x'\\ and 
the group operation of addition mod 1. Assume that the two-body potentials 
J(d(j,j')) and V(x, x'), G Z 2 , x,x' G S l , are of the form 




1, \j-j'\ = l, 

0, \j-f\^h 

— cos 2n(x — x'), p(x,x')<8, 
+oo, p(x, x') > 9, 



(3.2.1) 



with a usul agreement ■ (+oo) = 0, where 9 G (0, 1/4) is a constant. In this 
case, Hamiltonian H\ is equipped with Dirichlet's boundary conditions on 
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D C M A where 

D = |x A G M A : \x(j) - x(j')\ > 9 for some j,f G A with d(j,f) = 1 j. 

Then, V f3 G (0, oo), there exists an FK-DLR measure Jl — ftp G which is 
not 5 1 -invariant. Consequently, the corresponding FK-DLR state ip = (p^ G 
C5(/3) is not S* 1 -invariant. < 

Similarly to a ground state </?, we can define /J, a ground-state FK- 
measure. Namely, take an FK-DLR measure pL n G 0(/3 n ) an d consider its 
image ju n under projection u; r G 4 Xp 6 M r where x r = x r (u; r ) is 
the collection of initial points for loop configuration cj-p. Suppose that 
is a limiting point for sequence p, n as (3 n — > oo. Then we say that is a 
ground-state FK-measure. Furthermore, such a measure is called G-invariant 
if V finite set A C T, g G G and a (bounded) function : M A ° ->■ R, the 
integral ju(C/ A o(g)0) = £(0). 

Corollary 3.5. Suppose that the family of non-S 1 -invariant FK-DLR 
measures jlp specified in Theorem 3.4 has a limiting point ip as (3 — >■ oo. 
Then ^ is a non-S 1 -invariant ground-state FK-measure. < 
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4. Proof of the main results 



In this section we deliver proofs of the stated results. 

4.1. Proof of Theorem 3.1. Given finite sets A and A, A C A C T, 
the RDMK Ff (see (2.2.1)) is a continuous function on M A ° x M A °. The 
first observation is that the sequence Fw° n , is compact in C°(M A ° x M A °) 
by the Ascoli-Arzela theorem, as these functions are uniformly bounded and 
continuous. The latter property is based upon conditions (1.3.2)-(1.3.4) and 
the assumption that M is compact. 

More precisely, to show uniform boundedness, note that V finite A' C T 
with A' D A , the conditional RDMF g A °l r \ A ' (cf. Eqn (3.1.7)) satisfies 

g A °' r \ A '(aJV rvv ) < (e 2 ^)« A ° (4.1.1) 

for all path configurations uJ° G W^o, <jl>t\a> G Wr\A'- Here, the constant J 
is given by 



J = sup 



j'er 



(4.1.2) 



where J is the function from (1.3.3). (In fact, J coincides with the quantity 
J(l) in (1.3.1).) A similar bound holds if we replace T with ball A(n) C T\A'. 
After integration, this yields the estimate 

<)(x ,y )<(e 2 ^/)« A ° (4.1.3) 

where 

f = sup[pP(x,y),\V x pP(x,y)\,\V y pP(x,y)\ : x,y G M] . (4.1.4) 

Here p@(x, y) stands for the transition probability density from x to y for the 
Brownian motion (with the generator —A/2 where A is the Laplace operator 
on M) in time (3: 

^y) = -^ffm ^ exp (-|x-y + n| 2 /2/3), (4.1.5) 

n=(ni,-,n d )ez d 



The argument for uniform continuity (or equi-continuity) of RDMKs 

A(n)( 



F A ,V(x°,y°) is more technical. We want to check that the gradients 
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V x oFw\(x°,y°) and V y oFw\(x°, y°) are uniformly bounded. There are 
two contributions into the gradient: one comes from varying the measure 
P x o y o(daJ°), the other from varying the functional 

exp [— h A (ZJ° | u>A(n)\A )]- The first contribution can be uniformly bounded 
in terms of the constant p 13 . 

The second contribution can be analysed by deforming a path JUj G uf, 
j G A : one of the end-points x(j) or y(j) can be moved, say, along a geodesic 
on M (i.e., a straight line). The points on the path are then moved, at a 
scaled distance, via a parallel transfer (the afline connection on the torus is 
trivial). This contribution is related to the differentiation of the exponent 
and controlled due to the bound (1.3.2) on the derivatives of the potential. 

The second contribution yields an expression of the form 

/_ P x o, y o(du>°) ^2hj(u{j),u Hn) \ A o) 

JW *»,y° jeA o (4.1.6) 

xexp J-/i A °(aJ° | w A ( n )\A<0 



where functional hj(u>°, uJ\( n )\Ao) is uniformly bounded. Combining this with 
the argument used to estimate the RDMK F^(x°,y°) allows one to bound 
the gradients V x oF A ° n) (x°, y°) and V y oF A ° n) (x°, y°) as well. 

More precisely, given a site j G A , we need to differentiate in x(J) or 
y(j) the expression 



^2 /_ P x o, y o +n (rfa;)exp[-/;, AO (a;|a; A( „ ) \Ao] 
= ^ exp (-|x -y -n| 2 /2/3) 

n 

x / P x o iX o (dZJ°) exp [ - h A ° (p° + C„| w A („)\a°) 

Here we sum over vectors n = (n(j),j G A ) G (Z d ) A °. Furthermore, Cn is a 
linear map: 

Cn(r) = ^(y° + n-x ), 0<r</3. 

Finally, the measures P x o iy o +n and P x o ;X o refer to the standard Brownian 
motion on R d . 
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Suppose we differentiate in y(j) E R d . Differentiating the exponent yields 
a convergent series. Next, 



V y(j) exp 




Cn|^A(n)\A°^ 



T 



dT ^y(J) V )( u U, t),u(1, t)) 



x exp 



- h A (uf + Cn|^A(n)\Ao) 



is bounded due to (1.3.2). (The expression in the big brackets gives the term 
hj(u(j), ijt>A(n)\A°) figuring in Eqn (4.1.6).) This yields the desired result. 

Differentiation in x(j) can be done in a similar manner, by exchanging 
x° and y° in the above series. 



, A <> 



Now let an RDMK F be a limiting point for F A ,* 
Then we have that 



m 



C°(M A xM 



rA ^ 



lim 

k— >oo 



dx°dy° 



F A(n fc ) lX 



y°)-F A0 (x ,y ) 



n 2 



0. 



M A xM A 



A • uA° 

m ft 



In other words, the RDM R., 

' A(n fe ) 

RDM R A ° in the Hilbert-Schmidt norm: 
to Lemma 1 from 
well 



pA° pA° 



converges to the infinite-volume 
— > 0. According 



pA° pA° 



HS 



23] . the convergence takes place in the trace-norm as 
— > 0. We obtain that the sequence of states y?A(n) is 

tr 

w*-compact. 

In parallel, an argument can be developed that the measures /ia form a 
compact family as A /* T. More precisely, we would like to show that V 
finite A C T, the family of measures [i A ° is compact. To see this, it suffices 

L + l,L + 2,...} is 



to check that, for a fixed A , the sequence {/^[nj- 
tight and apply the Prokhorov theorem. 

To check tightness, we use the two facts, (i) The reference measure 
d^A° on W A ° (see Eqn (2.3.1)) is supported by loop configurations with 
the standard continuity modulus A/2eln (1/e). (ii) The probability density 
ao, c, _ d/$g(n(A(L))) 



n 



(cf. Eqn (2.3.3)) is bounded from above by a 



dp A{n) (Q(A(L))) 

constant exp {p [%A(L)]VT} (and from below by exp { - p [fJA(L)] V J*}). 
See (1.3.2) and (1.3.4). 
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The constructed family of limit-point measures /i A ' has the compatibility 
property and therefore satisfies the assumptions of the Kolmogorov theorem. 
The result is that there exists a unique probability measure \x on W T such that 
the restriction of /x on the sigma-algebra of subsets localized in A coincides 
with /i A °. 

The fact that /i is FK-DLR follows from the above construction. Hence, 
each limit point <p falls in class &{/3). □ 

Remark 4.1. Anticipating a forthcoming result for a general compact 
manifold M, we propose to discuss a version of the above argument for an ex- 
ample where M is a two-dimensional Klein bottle with a flat Riemannian met- 
ric. A convenient representation is through the universal simply connected 
cover which in this case is the Euclidean plane R 2 with the standard met- 
ric. For the fundamental polygon we take a square [—1/2, 1/2] x [—1/2, 1/2] 
where the following pairs of points (xi; x 2 ) are glued: 

(xi;-l/2) and (-xi; 1/2), where - 1/2 < x x < 1/2 

and 

(-l/2;x 2 ) and (l/2;x 2 ), where - 1/2 < x 2 < 1/2. 
The cover map T : R 2 — > M is as follows: 



T( Xl ,x 2 ) ^ ((-l) n2 ( Xl - m);x 2 - n 2 ) 
whenever - 1/2 < Xj - n x < 1/2, -1/2 < x 2 - n 2 < 1/2. 



(4.1.7) 



Here and below, n 1; n 2 G Z are integers. (In this example, G may be a circle 
S 1 (a ID torus) realized as an interval [—1/2,1/2] with points —1/2 and 
1/2 glued together. The action is: gx = (x x + g, x 2 ), for x = (xi;x 2 ), with 
addition in [—1/2, 1/2]. (Other choices of G (a ID torus of length 2 or a 2D 
torus) are analysed in a similar fashion.) 
In this example, the integral 

/ P x o iy o(dZJ°)exp [-/i A VVA(n)\A»)] (4.1.8) 



contributing to F j A °(x°,y°) can again be differentiated explicitly. For defi- 
niteness, take A to be a one-point set, say {o}, with particle configurations 
x° and y° reduced to single points in M (or rather in the fundamental poly- 
gon): 



x° = x = ( Xl ; x 2 ), y u = y = (y i; y 2 ), -1/2 < x h Yi < 1/2, i = 1, 2. 
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27T/3 



Then the above integral takes the form 

^ £ exp [-((-l)-x 1 -y 1 -n 1 ) 2 /(2^)] 

ez 

xexp [-(x 2 -y 2 -n 2 ) 2 /(2/3)] | Pf;f (dw? x du; 2 °) (4.1.9) 
-/i a ° (T[( Wl ° + 5; ni) ,u; 2 + ^ n2)N 



x exp 



k>A(n)\A° 



Here W^ 2 and Pf'^ 2 stand, respectively, for the space of continuous loops 
(closed trajectories beginning and ending at x) and the Brownian bridge 
distribution in the plane R 2 of the time-length (3. Next, 5^ are functions 
[0,(3] —> R providing deformations of plane loops w° from W£.'\. into plane 

paths Wl ° + 5i ni) fromWf: R y ' +ni : 



S[ ni \r) = T - (yi - Xl + m) , ^(r) = j (Y2 - x 2 + na) . (4.1.10) 

A similar formula holds after replacing W^ 2 and pf ^ 2 with W^y 2 and 
P^ 2 , mutatis mutandis. Differentiations V x and Vj, then become straight- 
forward (although rather tedious), confirming the above claim about uniform 
continuity of RDMKs (cf. Eqn (4.1.6)). 

Theorem 1.1 can be deduced from Theorem 3.1 since every limiting Gibbs 
state ip G (5° lies in (2» (in other words, 0° C 0), as follows from the above 
argument and the observations made in Sections 2.3 and 3.1. 

4.2. Proof of Theorem 3.2. We follow the approach initiated in [4]. 
The proof of Theorem 3.2 is based on the following bound for infinite- volume 
RDMFs: V finite A C T, uJ° e W A and g G G, 

q A \guf) + q A \g-'uf)>2 q A \uf). (4.2.1) 

Lemma 1 from [4] implies that (3.1.1) follows from (4.2.1). 

In turn, the bound (4.2.1) follows from a similar inequality for the con- 
ditional RDMFs: V finite A C T, uJ° G W A , g G G and a G (1, oo), for any n 
large enough and w r \A(n), 
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In fact, to deduce (4.2.1) from (4.2.2), it is enough to integrate (4.2.2) in 



d/z 



r\A(n) 



(a>r\A(n)) and let a \ 1. 



Now, Eqn (4.2.2) is deduced after performing a special construction re- 
lated to a family of 'gauge' actions g A ( n )\ A o on loop configurations u>A(n)\A°] 
see Eqns (4.2.4), (4.2.5) below. A particular feature of action g A ( n )\ A o is that 
it 'decays' to e, the unit element of G (which generates a 'trivial' identity ac- 
tion), when we move from A towards T \ A(n). Formally, (4.2.2) will follow 
from the inequality: V finite A C T, TU° G W\, g G G and a G (1, oo), for 
any n large enough, u> A ( n )\ A o and w r \A(n), 



■ exp 



a 

+-exp 



> 



h A{n) ((guJ°) V (g A(n) \ A oU> A(n) \ A o)|u> r \A(n); 

^((g" 1 ^ ) V (gX(„)\AO w A(n)\Ao)|wr\A(n)) 
" h A{n \cJ° V W A („)\Ao|Wr\A(n)) 



(4.2.3) 



exp 



Indeed, (4.2.2) follows from (4.2.3) by integrating in dz/ A ( n )\ A o(a; A ( ra )\ A o) and 
normalizing by £ A ( n \\ a° ( w r\A(n) ) ] cf. Eqn (3.1.16) with A' = A(n). (Here 
one uses the fact that the Jacobian of the map u>A(n)\A° l_ > gA(n)\A° a -'A(n)\A () 
equals 1.) 

Thus, our aim becomes to prove (4.2.3). The gauge family g A ( n )\ A o is 
composed by individual actions g^ G G: 



gA(n)\ A o = {gS* , jGA(n)\A°}. 



(4.2.4) 



Let us identify the element g G G with a vector 9 = OA G M and use the 
additive notation: gx := x + 9, x G M. Then G G corresponds to 
multiples of thye vector 9. Namely, we fix a positive integer value f such 
that A C A- and identify 



g^ with 9v(n,j) 

where 

1, d(o,j)<f 

^( d (j» -f,n-r), d(o,j)>f 
In turn, the function i?(a, b) satisfies 

l(0<a<6) ,! 



v(n,j) 



(4.2.5) 



(4.2.6) 



i?(a,6) = l(a < 0) + 



/ 2;(u)du, a, 6 G R, (4.2.7) 

«/ a 
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with the same functions Q(b) and z(u) as proposed in [4] 
Q{b) = [ z(u)du, 



(4.2.8) 

where z(u) = l(u < 2) + l(u > 2) — , b > 0. 

it In w 

Moreover, g^/ n s\ A o is the collection of the inverse elements: 

SA(n)\A« = {g? rl > 3 e A(n) \ A } . 

We will use the formulas for g^ for j G A(n), or even for j G T, as they 

agreed with the requirement that g^ = g when j G A and g^ = e for 

j G T \ A(n). Accordingly, we will use the notation g A(n) = 

Next, we use the invariance property (1.3.5). The Taylor formula for 
function V G C 2 yields for G A(n): 



V (gf c. 



(n) 



+r (gf ^g^ V) - 2y(w i ,o; i 
< C|0| 2 |<u(n,j) -v(n,f)\ 2 V. 



(4.2.9) 



Here C G (0, oo) is a constant V is taken from (1.3.2) and notations from 
(4.2.5) are used. 

The bound (4.2.9) is crucial and exploits the structure of the group action. 
It uses the fact that the first-order terms in the expansion in the left-hand 
side of (4.2.9) cancel each other because of the presence of elements gj and 

g^ and their inverses, g^ and . (This idea goes back to [14] and 

mo 
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The term \v(n,j) — v(n, j')\ 2 can be specified as 



\v(n,j) - v(n,j 



/\|2 



0, if d(j,o),d(f,o)<f, 
0, if d(j,o),d(j',o)>n, 
[#(d(j, o) — r,n — r) 
-tf(d(/,o) -r,n-r)] 2 , 

if f<d(j,o),d(f,o)<n, (4.2.10) 
^( d (j» -r,n-f) 2 , 

if r < d(j, o) < n, d(j', o) e]r, n[, 
tf(d(/,o)-r,n-r) 2 , 

if r < d(j', o) < n, d(j, o) e]r, n[. 

By using convexity of the function exp and Eqn (4.2.9), V a > 1, 



exp 



a 

+-exp 



> a exp 



- hHn) (gA(n) V W A(n)\Ao) |Wr\A(n)) 

^ W (gAW^ VW A W \Ao)|^A W ) 

- ^ A(n) (g A (n) (W° V WA(n)\Ao), W r \A(n) 



(4.2.11) 



-^ A(n) (gA(n) (55° V « A (n)VAo) l«r\A(»)) 



> a exp 



- /i A ^ (a; V u> A(n )\ A o | w r \A(nj 



-C*/2 



where 



* = *(n,g) = |^| 2 £ Ad(jj')) \v(nj)-v(njt 

(jj')6A(n)xr 

The next remark is that 



(4.2.12) 



*<3|#| 2 £ l(d(j,o)<d(j",o))j^ 

(j,j')eA(n)xr 



(4.2.13) 



x 



$(d(j, o) — r, n — r) — $(d(f, o) — r, n — r) 

where, with the help of the triangle inequality, for all j,j' : d(j, o) < d(f, o) 
< "&(d(j, o) —r,n — r) — $(d(y, o) —f,n — f) 



z(d(j,o)-r) 
~ U,J) Q{n-r) ■ 



(4.2.14) 
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This yields 



# < 



31101 



Q(n — r 



Y- \ 2 



J2 J(d(j,j'))d(j,f) 2 «j,o)-r) 2 



(jj')eA(n)xr 



W r; l^.,^ J , G A n+ro 

In view of (1.3.4) it remains to estimate the sum E z(d(j, o) — r) 2 . To 

jeA n+r() 

this end, observe that uz{u) < 1 when u G (3, oo). The next remark is that 
the number of sites in the sphere £„ grows linearly with n. Consequently, 

E z(d(j,o) - r) 2 = z(k-r) £ z(k-f) 

jeA(n+r ) l<k<n+r iGS fe 

<C J2 z{k - r) < dQ(n + r - r) 

l<k<n+ro 



and 



C 



— > oo, as n — > oo. 



Q(n - r) 

Therefore, given a > 1 for n large enough, the term 

ae -C*/2 in the RH g 

of (4.2.11) becomes > I. Hence, 



a 

2 6XP 
a 



+ -exp 



(§A(n) (W° V W A(n)\Ao) |w r \A(n)) 
hHn) (gX(n) ("° V «A(n)VA«) K\A(n)) 



(4.2.15) 



> exp - /i A ( ra ) (uJ° V w A(n)Vv o|u>r\A(n) 
Eqn (4.2.15) implies that the conditional RDMF 

? A ° |r ^ A(n) (^l w r\A(n)) = / dl/A(n)\Ao(w A (n)\Ao) 



X 



A(n)\AO 

exp [ - /j. A °(aJ° V ^A(n)\Ao|^r\Ao)] 

5A(n)(WT\A(n)) 



(4.2.16) 



obeys 



lim 

n—>oo 



>2 lim g A °l r \ A (")(ZJVr\ A(n) ) 

n— s-oo 

32 



(4.2.17) 



uniformly in boundary condition u>r\A(n)- Integrating (4.2.17) in 
d/iP A(n) (o;r\A(n)) yields (4.2.3). □ 

4.3. Proof of Theorem 3.4 and Lemma 1.1. In Theorem 3.4 our 
argument follows the idea proposed in [?]. On the lattice Z 2 consider the 
squares 

A(n) = {j = (ji,j 2 ) : max [|ji|, \j 2 \] < n}, n = 1, 2, . . . 

The outer boundary of A(n) is the set 

E(n + 1) = {j = (ji,j 2 ) ■ max \j 2 \] = n + 1}. 

Fix a point x* G 5 1 , a value (3 G (0, oo) and consider a state y?* = (p(jJ*) 
induced by measure /i* = /i^*) G (S(/3) which is a limiting point for the 
family of measures /x* = l^p^n) (2-3.10)) as n — > oo. Here /i* stands for 
the probability distribution on Wa(ti) with the 'cooled' boundary condition 

where 

u*{t) =x*, 0<t<(3. (4.3.2) 

Without loss of generality, assume that /i* = lim /x*. To simplify the 

notation, let us also omit the subscript *, writing // = //*. If state <p = <p* 
is not S' 1 -invariant, we are done. Otherwise, suppose that if is ^-invariant. 
Then choose an arc a = (x* — 1/200, x* + 1/200) of length 1/100 around the 
point x* and let n (a) be the orthoprojection on the subspace in "H — H 
formed by functions supported by arc a. Then 

^(Ilo(a)) = / d//{°>(a;o)l(xo(a;o) G a) = -L. (4.3.3) 
Jw {0} iUU 

(The lower/upper scripts and {0} indicate that we take A = {0}, i.e., 
consider spins attached to lattice site G Z 2 .) 

Hence, for n large enough, the conditional distribution 
d//{°}l E ( n+1 ))(a>o|ww n+1 )) for u G W{ }, given boundary condition wL n+I j, 
satisfies: 

f d/ i(WI 2 (" +1 ))(u;o|a;£ (n+1) )l(xo(u;o) € a) < L (4.3.4) 
Jw {0} Jy 
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Next, given rj G (0, 1], consider a family of points 

x jtV = x* +j 1 r]6 mod 1, j = (j u j 2 ) G Z 2 , 

and the family of the corresponding cooled loops {ojj^}'- 

u^{t) = x j>ri , 0<t<(3. (4.3.5) 

Further, consider the loop configuration uJj](n+i), v = j £ S(n+1)} over 

the boundary S(n+1) formed by loops uij tV - For r] — 1, the only configuration 
over A(n) compatible with the boundary condition u)Y,(n+i),ri is the one where 
all loops coincide with co v j: for any other choice of the configuration, the 
energy /i A ( n )l s ( n+1 ) is equal to +oo. By continuity, for each n, there exists 
rj(n) G (0, 1) such that the probability measures d/i^ }' s (™ +1 ^(w |^£(n+i),^(Yi)) 
conditional on £>£ n+1 ,rj(n)) satisfy 

f d/U (Wl s (" +1 ))(o;(0)|cD s(ri+1) ^ n) )l(a; ( W o) G a) = \, (4.3.6) 
Jw {0} 6 

and rj{n) is uniformly separated from and 1. Any limiting point Jl of the 
sequence of conditional measures /j( A ( n )l s ( n+1 ))( . |w S ( n+ i)^( n )), n — > oo, yields 

/ d^( Wo )l(x o mGa)4, (4.3.7) 

and the induced state (p gives ^(n a (0)) = 2/3. It means that neither fx nor 
are ^-invariant. □ 
Proof of Corollary 3.5. Eqn (4.3.7) guarantees that V (5 > there exists a 
non-S^-invariant measure fx^ G (25 (/3) . Passing to a limiting point as (3 — >■ oo 
yields a ground-state measure ^ with the property that 

/ d^( Wo )l(x Wea)4, (4.3.8) 
again contradicting S' 1 -invariance. □ 

Proof of Lemma 1.1. Let Ai > A2 > • • • be the sequence of the eigenvalues 
of operator R and ej(x), i = 1,2,... be the corresponding eigenvectors. As 
follows from (1.4.11), 



lim ^2 (( R nei, ej) - XiSij) = 0. 



1,3 
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We want to show that the sequence {R n } converges to R in the Banach 
space £ of the trace-class operators in L 2 (M, v) with the norm || ■ l^. We are 
going to use a natural basis in € formed by the system of rank one 'matrix 
units' Eij = |ej)(ej|. Set: 



p(*o) 
tin 



Next, set 



1<M<*0 i,j>k) 

Rn — ^2 ^2 (Rn&ii&j)Ejij. 
l<i<io j>io 



l<i<io 



Clearly, R^ and R^f are positive-definite operators. Furthermore, 



\R { J o) \ 



|if o) 



■n 1 1 1 



and 



R n = R^ + R { : o) +R^+(R^y . 
Take an arbitrary e > and choose i = i (e) and n = n (e) such that 



i>io 



and for n > 



n 



Then for n > n , 



\R-Rnh < || J R-i? (io) || 1 + R {io) -R { n o) 



p(*o) _ p 



< e/8 + e/8 + 
It remains to estimate the term 



p _ p(«o) 



< 



R 



<<o) 



p _ p(*o) 
p _ p(*o) 

-TL n -TL n 

= 1 - 



. To this end we write: 



(«o) 



< 1 - H^ )^ + e/8 + 2 
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p(*o) 
tin 



R 



(io) 



R 



(io) 



< 1 - e/4 + 2 



p(*o) 



Finally, 









1/2 




1/2 


o(*o) 


l<i<io 


^ ^(Rn&ii Gj) 
J>k) 


< «0 


^ ^{RnGj, Gj) 
- i¥=j 





which is < e/2. This completes the proof of Lemma 1.1. □ 
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